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Abstract 

In this paper, we present an automated system, whicii has the capabihty to catch 
and track solar hmb prominences based on observations from EUV 304 A passband. 
The characteristic parameters and their evolution, including height, position angle, 
area, length and brightness, are obtained without manual interventions. By apply- 
ing the system to the STEREO-B/SECCHI/EUVI 304 A data during 2007 April - 
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2009 October, we obtain a total of 9477 well-tracked prominences and a catalog of 
these events available online at http://space.ustc.edu.cn/dreams/slipcat/. A 
detailed analysis of these prominences suggests that the system has a rather good 
performance. We have obtained several interesting statistical results based on the 
catalog. Most prominences appear below the latitude of 60 degrees and at the height 
of about 26 Mm above the solar surface. Most of them are quite stable during the 
period they are tracked. Nevertheless, some prominences have an upward speed of 
more than 100 km/s, and some others show significant downward and/or azimuthal 
speeds. There are strong correlations among the brightness, area and height. The 
expansion of a prominence is probably one major cause of its fading during the rising 
or erupting process. 



1 Introduction 



Solar prominences, also called filaments when they are viewed on-disk, are long-observed 
but still not well-known structures in the solar atmosphere. Since they are outstanding 
features in multiple-wavelength observations of the Sun and have close relationships with 
various solar eruptive phenomena, prominences are always one of major topics in solar and 
space physics. The key issues in prominence/filament studies are their formation, main- 
tenance, dynamic processes, and their roles in other related solar activities, e.g., coronal 
mass ejections (CMEs) and flares. With t h e aid of modern sensor technology, man y facts 
of prominences are revealed [e.g., Poland . 19861 : Tandberg-Hanssen . 1995 : Martin . 1998 ; 
Patsourakos and Viaa . l2002l . and the references therein] . Prominences are dense (electron 
density ~ 10^ — 10^^ c m~^) and cool 5000 — 8000 K) plasmas floati i ig in hot and di- 



luted solar corona [e.g ., Engvold and BrynildseA . 19861 : Hiei et al . 1986; HirayaniA . 1986 



Madjarska et all Il999| . They can appear anywhere from active regions to polar regions, 
and live for days to months. They have spines and barbs, and always straddle above po- 
larity inversion lines. There are sometimes strong counterstreamings along spines and very 
dynamic vertical flows. The chirality of prominences/filaments obeys the pattern that most 
prominences on northern hemisphere are dextral while most ones on southern hemisphere 
are sinistr al. The assoc i ation rate of eruptive prominen ces with CMEs is more than about 



70% [e.g., \Gilbert et all . l2000l : iGopalswamv et al\ . 12003 



Many of the above findings are made through statistical investigations combined with 
case studies. A continuously updated catalog of prominences with unbiased parameters 
is undoubtedly helpful for such researches, especially in the age of the explosive growth 
of observational data. For instance, the successful launch of STEREO (Solar Terrestrial 
Relationship Observatory) spacecraft (A and B) in 2006 led the amount of solar observa- 
tions explosively growing to more than 12 GB a day, and now it has increased to about 
2 TB a day from SDO (Solar Dynamic Observatory), which was just launched in Febru- 
ary 2010. NOAA/SWPC0 routinely compiles a list of solar events, in which on-disk fil- 
aments and limb eruptive prominences are included; but the list is far from complete. 
Thanks to the unique properties of prominences/filaments, they can be clearly observed 
at multiple wavelengths, such as Ha, He I 10830 A, He II 304 A, radio waves, etc [e.g. 
Schmahl et al . 1974 : Hanaoka et oZI . Il994l : Penn et oZI . 19941 : Chiuderi Drago et al . 2001 



Labrosse and Gouttebroza . |2001| , and therefore it is possible to extract them from the vast 



amount of data automatically and consistently. 
^http:/ /www. swpc.noaa.gov/Data/index. html 
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Recognitions of on-disk filaments and limb promin ences are different. The former is 



mainl y accomplished by st udyi ng Ha data. For example, \Gao et al\ |2002| , \Shih and Kowalski 



2003j , \FuUer et al\ |2005j and \Zharkova et al\ |2005| developed codes to automatically de- 



tect fi laments in full-disk Hct images. The automated system developed hy lBernasconi et al. 
2005| is able to detect, classify and track Ha filaments efficiently. EUV observations are 
much more difficult to be used in detecting on-disk filaments due to the low contrast of 
filaments and the involvement of coronal holes. However, EUV observations are suitable 
for limb prominence detection. Through the usage o f Fe IX/X 171 A, Fe XII 1 95 A, Fe 
XV 2 84 A and He II 304 A images from SOHO/EIT \Delahoudmiere et al\ [l995|) instru- 



ment, \Foullon and Verwichtel |2006| developed algorithms to recognize limb prominences. 
In their method. He II 304 A data provide the basic criteria for the selection of candi- 
date prominence regions, and other emission lines are used to rernove ac tive regions, which 
also appear brightly in EUV 304 A. Most recently, iLabrosse et all 2010| also developed an 
efficient code to detect limb prominences in EUV 304 A images. 

In this paper, we present an automated system of detecting and tracking solar limb 
prominences based on only He II 304 A data, and a resultant on-line catalog as well, 
which can be continuously updated. The performance and limitations of the system are 
presented in section [Hi Based on our catalog, some preliminary statistical results of solar 
limb prominences are also presented. The reasons we choose He II 304 A emission line 
rather than Ha are the following. First, for prominence/filament observations. He II 304 A 
line is the only one uninterruptedly imaging the Sun with high cadence (operated by space- 
borne instruments SECCHI/EUVI on board STEREO twins, and AIA on board SDO). A 
complete database of limb prominences is therefore possible to be established. Secondly, 
the high time resolution of the data allows us to track their evolution, even small changes. 
Thirdly, the projection effect can be minimized for certain parameters, such as height, 
radial speed, etc. Fourth, they are complementary to the catalogs of on-disk filaments. 
Further, there is so far no well-established on-line catalog for limb prominences. 



2 Method 

Our system consists of five modules. The first module is to select prominence candidates; 
the second one is to extract necessary parameters for further usage; the third one is to 
discriminate prominences from other non-prominence features, such as active regions and 
noise; the fourth one is to track the prominences for the evolution; and the last one is to 
generate a catalog of prominences with final parameters. Here we use EUVI 304 A data 
from STEREO-B/SECCHI to illustrate these processes. 

2.1 Module 1: Prominence Candidate Selection 

The functions of module 1 are illustrated in Figure [H A raw EUVI 304 Aimage is shown in 
Figure [1^. The background brightness above the limb generally decreases with increasing 
distance, r, from the solar center. Similarly, the prominences near the solar surface are 
much brighter than those at high altitude (for example, comparing prominence A and B 
marked in the image). The variation of the brightness of the prominence with r is further 
discussed in Sec J4.3[ Although region B is too dark to be noticed in the raw image, we still 
consider it a prominence candidate for its higher density compared to the ambient coronal 
plasmas. 
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T he first step of proces sing is to use a tecfinique similar to tfie normalizing-radial-graded 
filter Morgan et al . l2006| to rescale the brightness so that the contrast is independent of 
r. To do this, a background image is first created, which is a circular symmetric image 
with respect to the center of the solar disk as shown in Figure [T|d. The pixel value at any 
r is just the average value of all the pixels along the circle at r in the original image. It is 
obvious that the brightness of background plasma does drop quickly as r increases. Then, 
we obtain the rescaled image (FigJTt) by using the following formula 



Rescaled Image 



Original Image (Fig{T^) + 6 
Background Image (Fig{T]D) + 6 



(1) 



Here, (5 is a small value to avoid dividing by a near-zero value. Both prominence material A 
and B become much clearer. For the STEREO-B/EUVI 304 A images, S is chosen to be 5 
through trial and error; however, this number may change for other instruments, depending 
on the signal-to- noise ratio. 



STERE0JB/EUVI304 




Figure 1: A sample image on 2007 October 8 to illustrate the processes of selecting promi- 
nence candidates, (a) The raw EUVI 304 A image, (b) circular symmetrical background 
image, (c) rescaled image, (d) binary image of selected kernels, (e) binary image of possi- 
ble prominence regions, and (f) the rescaled image with the boundaries of the recognized 
regions. The bright patch above the east limb in (f) is an active region, which will be 
removed by module 3. 

The further recognition, which applies the technique of region-growing with certain 
thres holds, is based o n the rescaled image. The fo llowing process is similar to those by. 



e.g., \Gao et al.\ |2002| and Wernasconi et al.\ |2005j . and thus we just briefly describe it 



here. First, we set a threshold thkni to pick all the pixels with larger value as kernels. The 
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searching region is from 1 Rs to Vmaxi, where Vmaxi is an upper boundary, above which 
there is no kernel selected. The boundary of Vmaxi is needed because the signal-to noise 
ratio will become low when approaching the edge of the telescope's field of view. For 
STEREO-B/SECCHI EUVI images, we choose the value of 1.7 Rs- The selected kernels 
serve as the seeds, from which the whole prominence regions grow out. Figure [T]i is a 
binary image showing the kernels. Here some small kernels, which are isolated pixels due 
to the presence of noise, have been removed by applying a morphological opening operator 
with a box size of Sn x Sn- Secondly, let these kernels grow by setting another smaller 
threshold t/ipro, i-e., all neighboring pixels with values larger than thpro will be included 
in the growing regions. For the cases that several regions close to each other but not 
connected, we use a morphological closing operator with a box size of SfYi X S^fi to merge 
them together. Some regions, whose areas are smaller than tharea, are discarded to further 
prevent noise-like features from being included. The resultant regions are the candidate 
prominences as shown in Figure [1^. Figure [T]F is obtained by superimposing the boundaries 
of the recognized regions on the rescaled image. The set of arguments discussed above are 
listed in Table [H 



2.2 Module 2: Parameter Extraction 

Once we have the boundary of a region of interest, the extraction of parameters of the 
region is straightforward. According to the scaling information in the header of FITS file 
of the image, we can calculate the area (A) and average brightness (F) of the region, the 
minimum and maximum positions in both radial and azimuthal directions {vhot, ^top, (^min, 
and 9max), and the centroid of brightness {rcen and Been)- It should be noted that, as the 
signal-to-noise ratio decreases significantly near the edge of the field of view, we set another 
upper boundary rmax2, slightly larger than Vmaxi, and consider that the parameters of any 
prominence extending into the region above it might be unreliable. 

Further, we linearize (i.e., get the spine of) the region to obtain certain morphological 
inform ation. Figure [2] p resents a sample. We use the morphological thin operator [refer 



to, e.g. \Lam et all . Il992| to get the skeleton (panel b) of the prominence of interest (panel 
a). Usually, a skeleton is too intricate because of many branches involved. To remove 
trivial branches, we first calculate the length (weighted by the rescaled brightness) of each 
branch. Then, for branches connecting to the same node, we compare their lengths and 
save the largest one. The above steps are iterated until there are only two ends (panel 
c). The resultant curve is further smoothed to get the spine (panel d) by applying the 2- 
dimensional mean filter method. One should note that any region of interest will be finally 
simplified to a line with only two ends even if it actually has three or more ends/footpoints. 

Since most prominences are loop-like structure in morphology, we keep the length of the 
spine as the characteristic length (L) of the recognized region. Meanwhile, the obtained 
spine can be used in 3D reconstruction of a prominence if it is viewed in two different visual 
angles at the same time (e.g., by the STEREO twins or combined with SOHO), which will 
be specifically studied in another paper. 

2.3 Module 3: Non-prominence Feature Removal 

While photons in EUV 304 A images ma inly come from He II emissions, they also have 



contaminations from hot coronal lines [e.s. JZhang et ami999| . As a result, prominences are 



not the only bright feature in EUV 304 A wavelength, and active regions also appear bright. 
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Figure 2: Extracting the spine of a prominence, (a) Original region, (b) Skeleton, (c) Spine, 
(d) Smoothed spine. 



6 



In \FouUon and Verwichtd |2006l | work, the authors reahzed this fact and used observations 
in other wavelengths to exclude active regions from their detected bright regions. The 
regions recognized through our first two modules also contain active regions and some 
noisy features. We do not, however, try to involve other observations in our detection, 
which will make the system more intricate and prone to additional errors. In our system, 
the previously extracted parameters for each recognized region will be used to discriminate 
real prominences from these non-prominence features, as discussed below. 

Prominences have a different appearance from other features. For example, in morphol- 
ogy, a prominence usually looks like a loop or stick, while an active region is shaped as a 
round blob. In brightness, a prominence is almost flat over radial distance in the rescaled 
image, while an active region is not. Thus one can use a classification method to remove the 
non-prominence features. There are many classification methods, e.g., linear discriminant 



analy sis (LDA), support vector machines (SVM), neural networks (NN) [e.g. AMeyer et al. 



2003| . The metho d we adopted here is the linear discriminant analysis. One can refer to 



the paper by, e.g., iFishen |1936l | for the principle of the linear discriminant analysis. 

Through many tests, the parameters \nA (standing for the size of a region), In^ (for 
the shape) and Inxp (for the variation in brightness, where Xf is the value of Chi-square 
goodness-of-fit statistic for the brightness F as a linear function of distance r) are chosen to 
construct the linear discriminant function (LDF). Our sample contains 5066 regions from 
a total of 3780 images (4 images per day, near 00:00, 06:00, 12:00, 18:00 UT, respectively, 
from 2007 April 1 to 2009 October 31). Each region is checked by eyes to determine 
which group it belongs to, the prominences or non-prominences. On the basis of this large 
collection of features of known classification, or the truth table, we derive the LDF as 

A 

X = 1.460 In A + 1.103 In 0.491 In (2) 

L 

or 



X = 14.20^^ + 4.70-|4v + 1-367!^ (3) 
{In A) (Inf) {InxD 



where the quantity (/), the mean value of / calculated based on our truth table, is used to 
normalize the parameters, and so that we can learn the importance of the parameters from 
their factors. The area A is the most important to discriminate a prominence from other 
features because it has the largest factor 14.20. However, it should be noted that some very 
big prominences might be missed due to the important role of area in the discrimination 
(for example the erupting prominence on 2009 November 2). But such missings usually 
take place in several frames, and therefore will not significantly affect the tracking result 
of the whole evolution process of the prominences identified in other frames. 

Figure [3] shows the discriminant result. It can be seen that the group of prominences 
(labeled Gl) generally have different LDF values from the group of non-prominence features 
(GO). Since there is still an overlap between the two groups, we evaluate the goodness of 
LDF by 

G = 1 - ^ (4) 
n 

where is the number of regions whose LDF value falls within the overlap, and n is the 
total number of regions in the truth table. In other words, the value of G is the ratio of 
the area of non-overlapped regions to the sum of the areas occupied by the two groups in 
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Figure 3: Result of the linear discriminant analysis of the truth table. The two groups, 
prominences (labeled Gl) and non-prominence features (GO), are indicated in red and blue 
colors, respectively. 

Figure [31 G = 1 means the LDF can completely discriminate the two groups. In our case, 
the goodness is about 0.86. 

Based on EqJH we can calculate the LDF value of any recognized region, and compare it 
with the derived distribution of the LDF values, which is fitted with Gaussian distribution 
functions as shown by the curves in Figure El to determine how likely the region is a 
prominence. The likelihood of a region being a prominence is given by 



mi 



mo + nil 



(5) 



where mo and mi are the values of the fitted Gaussian distribution functions corresponding 
to the LDF value for group and 1, respectively. A region with P < 50% is treated as a 
non-prominence feature and discarded. 



2.4 Module 4: Prominence Tracking 

Our method to track the evolution of a prominence is quite simple. Figure H] is the flow 
chart showing how we track a prominence. The top of the flow chart is a prominence to 
be tracked in an image, and the bottom of the chart gives the four possible results. Since 
the flow chart is detailed enough, we will not repeat it here. There are only several things 
that we would like to point out. 
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Figure 4: Flow chart to illustrate the prominence tracking process. 



First, the criterion used to judge if a prominence region is evolved from a prominence 
region in the previous image is to check whether or not there is an overlap between them 
in spatial domain. This requires that the cadence of the data should be high enough, 
especially when studying a fast erupting prominence. According to our statistical result 
(see Figj5]), which is plotted based on our catalog (refer to Sec J2.5p . most prominences move 
with a speed of about 4 km/s or less in radial or azimuthal direction, few ones may reach 
up to more than one hundred km/s. Considering that their characteristic length is L 60 
Mm, it is inferred that a cadence of 4 hours (for 4 km/s speed, or a cadence of 15 minutes 
for 100 km/s) is basically sufficient for prominence tracking, which is much longer than 10 
minutes, the cadence of STEREO/SECCHI EUVI 304 A data. 

Secondly, we use a time threshold thdis to determine whether or not a prominence has 
disappeared, i.e., if a previous named prominence has not been found in the successive 
images for a duration of thdis, it is treated as disappeared. A previous detected prominence 
may temporarily and intermittently 'disappear'. Such a 'disappearing' may not be real; 
it may be resulted from the unstable quality or jittering of images (although we have 
done some treatments on original data as mentioned in Sec J2.ip and/or small changes of 
prominence itself, which cause the brightness of the prominence to decrease down below the 
threshold th^ni or even thp^o temporarily. Note that this situation only happens to some 
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Figure 5: Histograms of the radial {left panel) and azimuthal speed {middle panel) of the 
centroid of prominences and the characteristic length [right panel). The average values are 
marked in the plots. The upper limits of the x-axes are chosen to make the plots readable, 
but do not mean the maximum values (The same treatment is made to the Figj6] and [7]) . 

small and/or faint prominences, not to major ones. Set a relatively long time duration thdis 
can efficiently track the entire evolution process of a prominence. Here we let thdis = 2 
hours. 

Thirdly, in our tracking process, the case that a prominence splits into two or more 
parts is considered (see the third result in the flow chart). However, we do not deal with 
the case of merging, which means there are more than one prominence regions (for example, 
A and B) in the previous image merging together and associated with only one prominence 
region (say C) in the current image. The merging of prominences is ambiguous, as the 
phenomenon can also be interpreted as that the prominence region A (or B) disappears 
and region B (or A) evolves to region C. In this scenario, no region merging takes place. 

Fourth, if a prominence is identified as a new one (see the left side of the fiow chart), we 
will check if it connects to the solar surface. Only those rooted on the Sun are considered 
as real prominences. This justification is based on the assumption that no newly emerged 
prominence is disconnected from the Sun. 

2.5 Module 5: Catalog Generating 

Solar prominences above the limb are identified and tracked by the above four modules. 
All the input arguments to SLIPCAT for the STEREO-B/SECCHI EUVI 304 A data are 
summarized in Table [H and the primary parameters extracted for each prominence at a 
certain time are summarized in Table HI Since we have the parameters of each prominence 
in its time sequence, it becomes possible to automatically extract its kinetic evolution 
information. For example, we can derive the velocity Vc and acceleration Hc of the centroid, 
from the variations of Vcen and 9cen with time. Also we can get some peak and average 
values of the above parameters, such as Amax, ^aue, etc. Moreover, for each prominence, 
we give the confidence level by the following formula 

r 1, P>_90% 
C = I 2, 75% < P < 90% (6) 
[ 3, 50% < P < 75% 

where P is the average value of the likelihood of the prominence over its period of tracking. 
A resultant online catalog is established at |http : /7 space . ustc . edu . cn/dreams/ slipcat/[ , 
where one can find all the final output parameters. The analyses in the following sections 
are based on the parameters in the catalog. 
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Table 1: List of the arguments used by SLIPCAT for STEREO-B/EUVI 304 A data 



Arp"nmprit^ 

1 U.lJ.lt'lJ.UO 




Units 


Tntpmrpt^i tinn 


6 " 


5.0 




A small value used in creating rescaled images (see 
Secl21j). 


fh, , 


9 




A friTPcinoln tot kpttipI ciplppfion i c;pp Spp 1 h 


fh °- 


1.7 




A threshold for region growing (see Secl2.1|). 


r 1 ^ 

' maxl 


1.7 




An upper boundary in r, above which there is no selected 
kernel (see SecJ2.1|). 


r ^ ^ 
' max! 


1 73 




An n Ti'M^^'r noi in H fi t^t in n t»si tsi in f^'f'f^'rc; r^t ti mr T^TOTni— 
u. |J iJCi Uwu.iiu.cii y 111 / . i lie [Jell dill" LCI o \JL cLii y ui vjiiii 

Tipripp P"v1"pri n 111 cr info fnp ypcrion i^noAT'P if ic; poncjinpypn fo 

IICIH^^C C-iVLCllLllll^ lllLw ullC 1 C^lWll CXUWVC lu lo L-'^llolU.Cl CJ. uW 

be unrehable (see Secj2.2|). 


S " 


5 


pixels 


Define a box used to remove noise-like kernels (see 
Sec 15^1). 


S c 


5 


pixels 


Define a box used to merge regions which are very close 
to each other (see SecJ2.1|). 


fh 

''area 


500 


Mm^ 


A threshold to remove very small regions (see SecJ2.1|). 




2 


hours 


A threshold to judge if a prominence has disappeared 
(see SecEl. 



Pixel values in rescaled images, probably changing for different instruments. 
^ Depending on the range of field of view and the signal-to-noise ratio. 
Need to be changed for different spatial resolution of images. 



Table 2: List of the primary parameters extracted for each prominence at a certain time 
Parameter Interpretation 

Coordinates of the centroid of the brightness. 
Give the span of a prominence in the radial direction. 
Give the span of a prominence in the angular direction. 
Area of a prominence in the units of Mm^. 
Characteristic length of a prominence in the units of Mm. 
Average brightness recorded in the original image. 
Likelihood of a recognized region to be a prominence. 



We note that it is straightforward to apply the system to other spacecraft data, e.g., 
SDO/AIA data (which is in our planning), though we present here only STEREO-B data. 
The only part we need to modify is the input arguments listed in Table [H Besides, the 
code is written in IDL language, and the modules 1-3 and modules 4-5 can run separately. 
The modules 1-3 do nothing with the causal relationship among the images, thus they can 
process images in serial or parallel. Using a computer with 2.33 GHz Inter Xeon CPUs, 3 
GB memory and Linux operating system, it takes 36 seconds on average for modules 1-3 
to process one full-size EUVI image and takes only 0.29 seconds for modules 4-5. Thus, 
to process the sequence of images over one day with the cadence of 10 minutes , it needs 
about 1.45 hours. The SDO/AIA data has higher resolution (4096 x 4096 pixels) and 
faster cadence (10 seconds), which will require modules 1-3 to spend much more time. It 
is estimated that it will take about 347 hours if we serially process such one-day images in 
the machine mentioned above. Thus, to reduce the processing time to less than 24 hours, 
there is a need of a small cluster with 15 or more CPUs to run the modules 1-3 in parallel. 



ycenj ^ceri) 
^hot ) T^top 
^minj ^max 

A 
L 
F 
P 
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which should be affordable. 



3 Performance and Limitations 

The SLIPCAT detects 19140 prominences from the STEREO-B data during the beginning 
of April 2007 to the end of October 2009. Figure E] shows the distributions of duration of 
the prominences and the number of frames in which a prominence is detected. It is found 
that 6348 prominences (33%) are detected in only one frame (upper panels), and the rest 
exhibits a roughly linear correlation between the duration and frame number as shown in 
the lower panel of FiglHl The solid line in the scatter-plot marks the expected relationship 
between the duration and frame number for the imaging cadence of 10 minutes. A point 
above the line means the instrument operates in a higher- cadence mode and the prominence 
is well tracked, while a point below the line implies that the prominence is missing in some 
frames. Most points distribute around the solid line. We define a prominence as poorly- 
tracked when it matches one of the following two criteria. (1) It is detected in only one 
frame, and (2) it is missing in 2/3 or more of the expected frames (marked by the dashed 
line in the scatter-plot). Note that the radial speeds of prominences are no more than 160 
km/s (see Fig|5]), which implies a prominence is expected in at least 5 or more frames. 
Thus a prominence detected in only 2 frames is also treated as poorly-tracked. There are a 
total of 9663 (50.5%) poorly-tracked prominences during the period of interest. As to the 
rest, we call them well-tracked prominences. 

These poorly-tracked prominences are generally small and their top portions (or leading 
edges) lie low. These are revealed in Figure [71 The average value of maximum areas of 
poorly-tracked prominences is 976 Mm^, about 3 times smaller than that of well-tracked 
ones; the average value of maximum lengths of poorly-tracked prominences is about 60 
Mm, nearly half of that of well-tracked ones; and the average top position of poorly- 
track ed prominenc es is about 21 Mm lower than that of well-tracked ones. As early as in 



1932. iPettitl |1932| had concluded that prominences are usually about 60 Mm long or more, 
10 Mm thick and 50 Mm high. These numbers are close to what we obtained here using 
modern data. By checking movies in the catalog, we find that the causes of such poorly- 
tracked prominences are probably resulted from the features to be marginal (close to the 
detection thresholds), contamination of nearby non-prominence features, and of course, 
the limitation of the detection algorithm. The parameters of these prominences may not 
be extracted correctly, and therefore they will be excluded from our statistical analysis in 
SecH 

The upper- left panel of Figure [6] suggests that some prominences may appear for up to 
260 hours, which seems too long to be possible. Due to the solar rotation, a prominence at 
a height of about 56 Mm (the average value of top position) above equator can stay visible 
for no more than 80 hours in theory. With the increasing latitude, this lasting duration may 
increase. If the prominence extends along the longitudinal direction, the lasting duration 
can be even longer. Figure [8] shows that the prominences with long duration generally 
appear in high latitude, where quiescent or polar crown prominences are usually present. 
Thus, it is possible to have such long-duration prominences. 

Figure |9] shows the daily counts of the well-tracked prominences. For most days, about 
14 prominences are detected. However, in extreme cases, there may be as many as 32 or 
as few as zero prominences a day. By checking the Ha images from BBSO, it is found that 
there are more filaments during solar maximum than solar minimum. During 2007 - 2009, 
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Figure 6: ( Upper panels) Histograms of the duration and the number of frames. The 
average values are marked in the plots. {Lower panel) Relationship between the duration 
and frames. The solid line indicates the expected relationship at a cadence of 10 minutes. 
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Figure 7: Histograms of the maximum area {left panel), characteristic length {middle panel) 
and top position {right panel) for the well-tracked (solid line) and poorly-tracked (dotted 
line) prominences. The average values (red for well-tracked and blue for poorly-tracked) 
are marked in the plots. 



the extreme solar minimum, there are usually only a few filaments in an Ha image, that is 
inconsistent with our results. To make sure that the prominences identified by SLIPCAT 
are not non-prominence features, we compare the EUVI 304 A images with the Ha image 
as shown in Figure [TOl The date of these images is 2009 October 7th, which is arbitrarily 
chosen; on this day, the STEREO-B spacecraft was 58 degrees behind the Earth. Thus 
the west limb in STEREO-B/EUVI 304 A image corresponds to about 33 degrees west to 
the central meridian in the Ha image. One can find that the western hemisphere in the 
Ha image is largely free of features except there are three prominences standing at high 
latitude (as denoted by arrows). Actually, there is a prominence at the low latitude, which 
can be clearly seen in the EUVI images (marked by a circle). The comparison demonstrates 
that there are probably some small prominences visible in EUV 304 A emission line, but 
invisible in Ha. It is also supported by many observati onal facts tha t solar filaments are 



generally more extended in EUV lines than in Ha [e.g., \Heinzel et all |2001| . Thus, it can 
be concluded that SLIPCAT is sensitive to recognize prominences, even those invisible in 
Ha . 

Further, the accuracy of the parameters listed in the catalog should be addressed. First 
of all, it should be noted that the values of these parameters just give us the information 
of prominences during the period they are detected, not necessarily through their whole 
life-time. Second, some parameters suffer from the projection effect, e.g., the value of 
area depending on the angle of view. Third, some parameters, e.g., velocity, change rates 
of area and brightness, are automatically derived from fitting data points. Thus their 
accuracy depends on the number of frames, the fitting function and the complexity of 
the prominence. As an example. Figure [TT] presents an erupting prominence observed on 
2008 November 13. This prominence rose from the south-east limb, and partially erupted. 
During its eruption, the prominence split into 3 major parts (see FigJTTli). one escaped 
from the Sun, one remained on the Sun, and the other erupted but returned back to the 
Sun later. SLIPCAT tracks the whole eruption process. The left panel of Figure [T2] displays 
the height-time profile of the leading edge of the prominence. The solid line is a quadratic 
fit through the data points and the dashed line is a linear fit. Due to the splitting of the 
prominence, the fitting curves obviously do not reflect the reality. The correct treatment is 
to study the evolution of the split parts separately. The right panel shows the height-time 
profile of the escaping part. However, the fitting result for that part is still not satisfactory 
enough. The time at about 11:00 UT is a critical point, before which the escaping part 
was slowly rising, while after which it experienced a fast acceleration and erupted quickly. 
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Figure 8: Histograms of the latitude for long (> 80 hours, dotted line) and short (< 80 
hours, solid line) duration prominences, respectively. The average values (red for short- 
duration and blue for long-duration) are marked in the plots. 
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Figure 9: Distribution of daily counts of limb prominences. 
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Figure 10: Comparison of the EUV 304 A prominences viewed by STEREO with Ha 
filaments viewed from the Earth. From the left to right, they are orig inal EUV 304 A image, 
rescaled image, rescaled image with selected prominences, and Ha image from BBSO. 

Thus a two-stage linear fitting is more appropriate than a pure linear or quadratic fitting. 
However, we still choose linear and quadratic functions to fit all the detected prominences, 
because this treatment can be easily operated in an automated way. One should keep it in 
mind that the fitting results only give a coarse estimation of the average speed (or change 
rates of area and brightness) and acceleration. 

4 Preliminary Statistical Results 

Now SLIPCAT has a complete dataset for STEREO-B data. A catalog of prominences 
seen from STEREO-A will be generated soon. Here we would like to present some statis- 
tical results of the STEREO-B prominences. Since we probably do not get very accurate 
parameters for prominences as discussed in the last section, the results obtained here are 
just preliminary. However, from a statistical point of view, these results should be sig- 
nificant. In the following analysis, all the poorly-tracked prominences are removed. Such 
prominences are generally extremely small and stay at low altitude as discussed in Secl3l 
thus the statistical results might suffer from the bias of selection; one should treat it as a 
statistics of moderate and major prominences. Moreover, we will include prominences with 
all confidence levels. There are only 257 (2.7%) prominences with confidence level of 2 or 
3, which represents a rather small fraction in the database . 

4.1 Static Parameters 

First of all, static parameters of the 9477 well-tracked prominences are investigated. These 
parameters are (1) the heliographic latitude of the centroid of a prominence at the first 
detection, (2) the mean value of the height of the centroid, (3) the mean value of the area 
and (4) the mean value of the brightness of a prominence. The distribution of heliographic 
latitudes (folded at equator) shown in the upper-left panel of Figure US] suggests that 99% 
prominences appear below 60 degrees. One may suspect that the lower count at high 
latitude may be caused by the perception that a prominence at high latitude is usually 
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Figure 11: An erupting prominence on 2008 November 13, which spht into three parts 
during the eruption. One escaped from the Sun, one remained on the Sun, and the other 
erupted but returned back to the Sun. 
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Figure 12: Height-time profiles of the leading edge of the prominence presented in FigJTTl 
The left panel is for the whole prominence system, while the right panel is for the escaping 
part only. The solid line is a quadratic fit through the data points and the dashed line is 
a linear fit. 
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Figure 13: Histograms of static parameters: [upper-left panel) heliographic latitude at the 
first appearance, mean values of [upper-right panel) height of centroid, [lower-left panel) 
area and [lower-right panel) brightness. The average values of these distribution are marked 
in the plots. 
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quiescent or a polar crown prominence and generally is long lived and extended. However, 
the scatter-plot of the latitude versus duration in the upper-left panel of Figure HH indicates 
that this perception is wrong, as it clearly shows that the durations of the prominences 
above 60 degrees are short. The long-duration prominences appear around 30 to 60 degrees, 
which probably implies that long extended prominences arise there. Similar result can also 
be found in Figure [H] though there we include the poorly-tracked prominences. Moreover, 
from rescaled EUV 304 A images, one can clearly find that the regions above 60 degree 
are generally occupied by polar coronal holes. By checking the catalog and movies, we find 
that the detected 'prominences' above 60 degrees are usually polar jets. Since their number 
is small, the statistical results obtained below will not be affected by including these 'false' 
prominences. 

The distribution of the heights of centroids suggests that about 82% of prominences stay 
at around 26 Mm above the solar surface. There is no obvious dependence of the height on 
the latitude as shown by the scatter-plot in Figure [HI The previous statistical study by 



Ananthakrishnan 


Mm. 


Pettit 1932 



196ll | suggested that the heights o f prominences var y between 15 and 150 



gave a value of about 50 Mm, and lKim et al\ |1988[ | showed that there is 



a peak at about 28 Mm in the distribution of heights. Theoretical work suggests that the 
height of a promi nence depends on t he gr adient of ambient coronal magnetic field strength 



with height [e.g., iFilippov and Peru . l2000l |. Our result is consistent with these studies. 
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Figure 14: Correlations of heliographic latitude with {upper-left panel) duration, mean 
values of {upper-right panel) the height of centroid, {lower-left panel) area, and {lower- 
right) brightness. 

The area and brightness can be usually used to evaluate if a prominence is a major one 
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or not. Similarly, we use the mean values of them to show the distribution. The average 
projected area on the plane of sky of all prominences is about 1072 Mm^. Nearly 60% 
prominences have smaller area than the average value. Figure [T3] suggests that the area 
is unrelated with the latitude. The brightness is recorded as digital number (DN) by the 
CCD camera. Its distribution is close to a Gaussian one, with an average value at around 
275 DN. It can be as low as 20 DN or as high as 960 DN. The scatter-plot in Figure [T^ 
shows a weak dependence of the brightness on the latitude. The low-latitude prominences 
may reach to a higher brightness than middle to high-latitude prominences. 



4.2 Dynamic Parameters 

For dynamic properties, we investigate the velocities (both radial and azimuthal) and 
change rates of area and brightness of prominences. Here we use the leading edge rather 
than the centroid in the analysis of the radial speed. One also may use the radial speed of 
the centroid, but it will bring a large error in the case that a prominence splits into two 
parts, one erupting and the other staying on the Sun. The histograms in Figure [T5] show 
the radial and azimuthal speeds. A speed of 10 arcsec/s corresponds to about 33 km/s 
at solar surface. More than 80% of prominences have no obvious motion in either radial 
or azimuthal direction, and the value of zero is the most probable speed (as shown in the 
insets). A few prominences may move upward at about more than 100 km/s, and also there 
are 37 (~ 0.4%) prominences having a radial speed < —20 km/s or an azimuthal speed 
> 10 arcsec/s. The former could be easily understood that an erupting prominence may 
have a large outward speed, but the cause is not obvious for the latter. 

We have mentioned in Secj3] that the quality or precision of the fitting results largely 
depends on the number of measurements. More measurements can efficiently reduce errors, 
and therefore the fitting results will be more reliable. By checking the movies of the 37 
prominences with the unusual speed, we find 34 prominences are detected in no more than 6 
frames, and the speeds of most (not all of) these 34 prominences are not correctly reflected 
by the fitting. However, since we have a large sample in the statistics, such a small fraction 
of corrupted events will not distort the overall picture shown in Figure [151 On the other 
hand, it is also realized that there are indeed some prominences having a large downward 
or azimuthal speed. The large downward speed may either present a real motion or is just 
resulted from the shrink of the prominence. Any further analyses on such extreme events 
will be pursued in the future. 

Similarly, for most prominences the change rates of area and brightness are quite small 
although the average value of the change rate of area is about 1.35 x 10^ Mm^/s. There 
are 106 (1.1%) prominences with absolute value of change rate of area > 10^ Mm^/s or 
brightness > 0.25 DN/s. The movies reveal that some prominences do change this fast. 



4.3 Fading of Prominences 

It is well known that prominences generally b ecome dimmer as they ri s e. The reason could 



be the heating of prominence r nateri als \e.s.. \Mouradian and Martresl. Il986 : 



19981 : iHanaoka and Shinkawm. 1199911 . mass loss [e.g., iRusin and Rybansky 



Ofman et all 


1982 


and/or 



expansion [e.g., \Bemporadl l2009l |. The first one belongs to thermal p rocesses, and the 
other two are dynamic proces ses [e.g., Mouradian and Martres . 1986; Mouradian et a/.l . 



19951 : 1 Tandherg-Hanssem Il995l | . Here we will look into this issue in a statistical way. In our 



parameters, we have the information of altitude, brightness and area of prominences, and 
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Figure 15: Histograms of dynamic parameters: {upper-left panel) radial speed of leading 
edge, {upper-right panel) azimuthal speed of centroid, {lower-left panel) change rate of area 
and {lower-right panel) brightness. The averaged values are marked in the plots. 

have no direct information about the heating or mass of prominences. Thus it is impossible 
to make a comprehensive study of the causes of the prominence fading, but we can learn 
how significant the factor of the expansion is. 

The first panel in Figure [16] shows a strong anti-correlation between the brightness 
and the height of prominences. The higher altitude makes the prominence dimmer. The 
diamond symbols with error bars indicate the average values of brightness around certain 
altitude. These points are fitted to obtain the following empirical formula 

F = 7.47 X 10^{h + 1.29)-°-^^^ DN (7) 

where h > 35 Mm is the height from solar surface in units of Mm. It describes the 
dependence of the brightness on the altitude. The points below 35 Mm are excluded in the 
fitting as they seem to follow another pattern. 

The second panel exhibits an evident positive linear correlation between the area and 
height. It means a prominence at a higher altitude tends to be larger. This phenomenon 
supports the picture that, when a prominence rises or erupts, it expands as well. The 
expansion is probably caused by the weaker constraint of ambient atmosphere at a higher 
altitude. Similarly, we fit the data points marked by the diamond symbols, and get 

A = 64:h- 906 Mm^ (8) 
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The reversed correlations of EqJT] and M suggest that the expansion of prominences must 
be a cause of the prominence fading when they are rising or erupting. 
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Figure 16: Scatter plots between minimum brightness, maximum area and maximum height 
of leading edge. 

Combining the two equations, we derive the relationship between the brightness and 
area (indicated by the dashed line in the last panel of Fig JT6|) 

F = 3.04 X 10\A + 9.88 x lo^)-^-^^! DN (9) 

where A is in the units of Mm^. Also we can directly fit the data in the last panel that 
leads to 

F = 3.96 X 10^{A + 2.42 x iQ^y^-^^^ DN (10) 

as shown by the solid line. Approximately, area is proportional to y^/^ where V is the 
volume of a prominence. If there is no mass loss or gain, the density of a prominence p 
is inversely proportional to the volume, i.e., p oc V"^ oc The brightness can be a 

proxy of the density of the plasma in the temperature window corresponding to EUV 304 A 
emission line. If there is no heating or cooling, the function F = co{A + ci)~^'^ is supposed 
to describe the relationship between the brightness and area. By fitting the data points, 
we obtain the dotted line given by 

F = 4.76 X 10^{A + 1.66 x lO^)"^'^ DN (11) 
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It is found that the sohd and dashed hnes are close to the dotted one, which imphes that, 
in a statistical point of view, the expansion is probably one of the major causes of the 
fading of prominences during their rise or eruption. Of course, t his statistical conclu sion 
could not be true for all individual cases. As revealed by, e.g., lOfman et al\ jl998| and 



Hanaoka and Shinkawd [1999|, the heating process may play an important role in the dis- 



appearances of some prominences. 



5 Summary 

We have developed an automated system of catching and tracking solar limb prominences in 
EUV 304 A images. The system, called SLIPCAT, is able to generate (1) a catalog of solar 
limb prominences and (2) some characteristic parameters of each detected prominence, in- 
cluding the height, position angle, area, length, brightness and their first and second deriva- 
tives with respect to time. SLIPCAT is composed by five modules, (1) prominence candi- 
date selection, (2) parameter extraction, (3) non-prominence feature removal, (4) promi- 
nence tracking, and (5) catalog generating. At present, an online catalog for STEREO- 
B/EUVI 304 A data has been generated (refer to|http : //space .ustc . edu.cn/dreams/slipcat/ ), 
and catalogs for STEREO- A/SECCHI/EU VI and SDO/AIA data are in preparation. 

Based on the STEREO-B/EUVI 304 A data, SLIPCAT proved to perform well in 
detecting limb prominences. 

1. It can distinguish real prominences from non-prominence features, e.g., active regions, 
without observations in other wavelengths by using the technique of linear discrim- 
inant analysis. The goodness of the overall classification is at about 86% successful 
rate. 

2. It detects as many as 9477 well-tracked prominences during 2007 April - 2009 Octo- 
ber, which means in average about 10 events per day. Compared to Ha data, it is 
found that SLIPCAT is sensitive enough to recognize almost all prominences, even 
those invisible in Ha images or very small ones. 

3. Thanks to the high-cadence EUV 304 A data, SLIPCAT is able to provide the detailed 
evolution processes of prominences quantitatively without manual interventions. The 
upper-right panel of Figure E] implies that a well-tracked prominence is at least de- 
tected in 28 images in average; and the case in Figure [12] shows that such high-cadence 
detection allows us to make a detailed analysis of its evolution, including its eruption, 
oscillation, etc. 

However, not all the parameters extracted by SLIPCAT can precisely reveal the real be- 
havior of prominences. The limitations have been addressed in the last paragraph of Secl3l 
Summarized here, they are that (1) the parameters characterize the properties of promi- 
nences during the period they are detected, not over the whole life-time, (2) they suffer 
from the projection effect, and (3) the speeds, change rates of area, length and brightness, 
which are derived from linear and quadratic fittings, may not be accurate. 

By applying SLIPCAT to the STEREO-B/EUVI 304 A data from 2007 April to 2009 
October, we obtain the following preliminary statistical results of solar limb prominences. 

1. In average, there are about 10 prominences standing above the solar limb per day 
during the solar minimum. For most days, about 14 prominences are expected to be 
detected, and sometimes the number could be as large as 32 or as small as zero. 
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2. Most (99%) prominences appear below latitude of 60 degrees, and the long extended 
prominences tend to arise between latitudes of 30 and 60 degrees. 

3. Most (82%) prominences have a height of about 26 Mm from the solar surface. 

4. The projected area of a prominence on the plane-of-sky is about 1072 Mm^ in average, 
and nearly 60% of prominences have a smaller area. 

5. Most prominences are quite stable during the period they are detected; no obvious 
change in position, area or brightness can be found. 

6. Particularly, more than 80% of prominences do not show obvious motion in either 
radial or azimuthal direction. However, some prominences have an upward speed 
of more than 100 km/s, and a few prominences present a significant downward or 
azimuthal speed. 

7. The brightness of prominences is anti-correlated with the height. The prominences 
at higher altitude look dimmer. 

8. The area of prominences is positively correlated with the height. The prominences 
at higher altitude are generally larger. 

9. From the statistical point of view, the expansion of prominences is probably one of 
the major causes of the fading of prominences during their rise or eruption. 
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